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ABSTRACT 

We present the results of a 10 ks XMM-Newton observation of NGC 3998, a "type-I" LINER galaxy 
(i.e., with significant broad Ha emission). Our goal is to test the extent to which radiativcly-inefficient 
accretion flow (RIAF) models and/or scaled-down AGN models are consistent with the observed 
properties of NGC 3998. A power-law fit to the XMM-Newton spectra results in a power-law slope 
of r = 1.9 and 2-10 keV flux of 1.1 x 10~^^ erg cm~^ s~^, in excellent agreement with previous hard 
X-ray observations. The OM UV flux at 2120A appears to be marginally resolved, with ~ 50% 
of the flux extended beyond 2". The nuclear component of the 2120A flux is consistent with an 
extrapolation of the X-ray power-law, although ^ 50% of the flux may be absorbed. The OM U flux 
lies significantly above the X-ray power-law extrapolation and contains a significant contribution from 
extragalactic emission. The upper-limit for narrow Fe-K emission derived from the XMM-Newton 
spectra is 33 eV (for Ax^ — 2.7). The upper-limit for narrow Fe-K emission derived from a combined 
fit of the XMM-Newton and BeppoSAX spectra is 25 eV, which is one of the strictest limits to date for 
any AGN. This significantly rules out Fe-K emission as is expected to be observed in typical Seyfert 
1 galaxies. The X-ray flux of NGC 3998 has not been observed to vary signiflcantly (at > 30% level) 
within the X-ray observations, and only between observations at a level of ~ 50%, which is also in 
contrast to typical Seyfert 1 galaxies. The lack of any reflection features suggests that any optically- 
thick, geometrically-thin accretion disk must be truncated, probably at a radius of order 100-300 (in 
Schwarzschild units). RIAF models fit the UV to X-ray spectral energy distribution of NGC 3998 
reasonably well. In these models the mid-IR flux also constrains the emission from any outer thin disk 
component that might be present. The UV to X-ray SED is also consistent with a Comptonized thin 
disk with a very low accretion rate {M < 10~^MEdd), in which case the lack of Fe-K emission may be 
due to an ionized accretion disk. Accretion models in general do not account for the observed radio 
flux of NGC 3998, and the radio flux may be due to a jet. Recent jet models may also be consistent 
with the nuclear fluxes of NGC 3998 in general, including the X-ray, optical/UV and mid-IR bands. 
The (ground-based) near-IR to optical photometric data for the nuclear region of NGC 3998 contains 
large contributions from extra-nuclear emission. We also derive nuclear fluxes using archival HST 
WFPC2 data, resulting in meaningful constraints to the nuclear SED of NGC 3998 in the optical 
band. We discuss a possible OM U band and USNO-B detection of the NGC 3998 ULX. 
Subject headings: X-rays: galaxies — galaxies: individual (NGC 3998) — accretion, accretion disks 



1. INTRODUCTION 

Since LINERs (low-ionization nuclear emission-line re- 
gions; Heckman 1980) were first identified as a class the 
nature of the source of ionization has been debated, with 
the two main possibilities being a starburst or AGN nu- 
clear source. If an AGN is responsible for the ioniza- 
tion, the next question is how similar are these AGN 
to Seyferts and QSOs. An interesting development has 
been the d iscovery of broad Ha emiss ion in many LIN- 
ERs po~F ihppcn ko. k Sarge'^IT997|) . These LINERs 
are almost certainly AGN powered, particularly when 



they are observed in early-type galaxies where outflows 
from massive stars and/or supernovae are not expected 
to be relevant. Early- type galaxies also tend to har- 
bor very massive central black holes, making broad- line 
LINERs observed in early-type galaxies ideal laborato- 
ries for studying (potentially) low accretion rate physics. 
Of course, this is assuming that the X-ray emission is 
isotropic and scales with accretion rate, however alterna- 
tively the accretion energy may be dissipated anisotropic 
and/or non-radiatively (i.e., in outflows; Di Matteo et al. 
2000). 
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NGC 3998 is an S O galaxy at a distance of 14.1 Mpc 
l)Tonrv et al. ll20?)]l) . Its Ha line contains a broad com- 
ponent with 37% of the Ha + [N II] flux and a FWHM of 
- 2000kms-i (Ho. Fih ppenko. fc Sargent 1997j. ASCA 
and iJeppoS'^X observations find that the 2-1 keV lumi- 
nosity of NGC 3998 is ~ 3 X 10^^ erg s'^ (.Pellegrini et all 
l2000HTerashima et al.l l2002'). and the 0.1-100 keV lumi- 
nosity is ~ 1 X 10^^ erg s~^. The mass of a putative 
black hole in the nucleus of NGC 3998 is estimated to 
be ~ 1O^M0, giving Lx/LEdd 1 x 10~^, or well into 
the radiatively-inefficient accretion flow (RIAF) regime. 
The main goal of this paper is to contrast NGC 3998 
with "normal" Seyfert 1 galaxies and to test the applica- 
bility of current RIAF models to the observed properties 
of NGC 3998. 

We present the results of a recent XMM-Newton ob- 
servation of NGC 3998, analyzed in conjunction with 
archival X-ray data. Previous X-ray observations showed 
that no Fe-K line is present in the X-ray spectrum from 
NGC 3998, and here we tighten this constraint, which 
places limits on the accretion geometry in NGC 3998. 
Additionally we will use XMM-Newton Optical Monitor 
data along with other multi-wavelength data to compare 
the spectral energy distribution (SED) of NGC 3998 with 
those predicted by various accretion models. 

2. DATA REDUCTION 

The parameters of X-ray observations are listed in 
Table [T] The ASC A d ata was redu ced as described 
in iPtak et all (|1999|) and iTerashima et all (|2002|) . The 
BeppoSAX data were reduced by the XSELECT and 
FTOOLs in HEASOFT 5.2 from the cleaned event 
files obtained from the BeppoSAX Science Data Cen- 
ter. LEGS and MEGS spectra of NGC 3998 were ex- 
tracted from a circular region with a 4 arcmin radius. 
Background spectra were taken from an annular region 
around NGC 3998. Here we present CCD spectra from 
the three EPIC detectors (two MOS and one PN) and the 
Optical Monitor U-band and UV (UVW2 filter) images. 
The XMM-Newton data was reduced using XMMSAS 
version 5.4.1 with the use of XAssist, a software pack- 
age t hat performs the initi al steps of X-ray data anal- 
ysis llPtak & GrifRthsll200l . The XMM PN and MOS 
spectra were extracted from 37" circular regions (derived 
from fitting a Gaussian surface brightness model to the 
X-ray images). We compared the radial profile of the 
mosl image (chosen since the target was not close to a 
ccd boundary in that case) to a SciSim simulation using 
the spectral model (see i)3.1|l and found that the X-ray 
emission from NGC 3998 is dominated by an unresolved 
source. The only other significant source in the vicin- 
ity of NGC 3998 is NGC 3998 X-1 IIRoberts fc Warwickl 
1200011 which is discussed in The OM data was re- 

processed with XMMSAS and the images and source lists 
produced by the tool "omichain" were used in this paper. 
The OM spatial resolution depends on photon energy 
and the FWHM of the point-spread function (PSF) is 
1.6" and 1.9" in the U and UVW2 bands {XMM-Newton 
User's Guide). All spectra were binned to 20 counts/bin 
to allow the use of the statistic. 

3. RESULTS 
3.1. Spectral Analysis 



We fit the XMM-Newton (CCD) spectra of NGC 3998 
with a simple absorbed power-law model, allowing the 
overall normalization of the model fit to each detector 
to be independent. This resulted in a good fit with 
F ^ 1.9, consistent with the previous ASCA and Bep- 
poSAX ohseivations as shown in Table|21 Note that while 
the three observations resulted in very close agreement 
for the power-law slope, the column density measure- 
ments varied from 3.3 x 10^° cm~^ (or a factor of ~ 3 
greater than the Galactic value; Murphy et al. 1996) 
to 8.0 x 10^° cm^^. Given the calibration uncertainties 
below 0.5-1.0 keV for these detectors it is not clear to 
what extent these variations are significant. We show the 
power-law fit to the XMM-Newton data in Figuren(the 
BeppoSAX and ASCA data are plotted in Pellegrini et al. 
2000 and Terashima et al. 2002). The ASCA data also 
resulted in a flux ~ 30% lower than the BeppoSAX and 
XMM-Newton measurements (see ii3.3|l . Note that the 
PDS flux was ~ 45% higher the mean of the other flux 
measurements, however the bandpass of the PDS data 
was ~ 10 — 100 keV so the 2-10 keV flux is based on an 
extrapolation. We also tried restricting the BeppoSAX 
data to only include data below 10 keV (i.e., excluding 
the PDS data) and found no significant impact on the 
fit parameters. Power-law fit parameters from a fit to 
the XMM-Newton and BeppoSAX data simultaneously, 
with the normalization of each data set allowed to vary, 
is shown in Table [3 and only increased by 3 when the 
normalization of the PDS component was fixed result in 
the same 2-10 keV flux as the MFCS flux (i.e., the PDS 
data systematically deviate from the power-law fit at a 
confidence of only ^ 90%). 

There are no systematic residuals suggesting additional 
spectral features might be present. We do not find evi- 
dence for any edges suggesting that an ionized absorber 
is not present (note that the edges discussed in Pelli- 
grini et al. were detected at less than the 3a level). The 
RGS data do not have high enough signal-to-noise to con- 
strain spectral features. There is a slight excess of counts 
in the 0.65-0.7 keV region but this feature is not statis- 
tically signifi cant. We a lso tried a plasma plus power- 
law fit (|Ptak et all 119990 to the XMM-Newton spectra 
and found a significant improvement in x'^ ix^ /dof = 
1406/1334 compared with 1423/1338 for the power-law 
fit). An absorption component applied to the whole spec- 
trum resulted in Nh = 4.6(1.0 - 5.8) x 10^° cm~^ and a 
separate absorbed applied to just the power-law compo- 
nent resulted in an upper limit of 4.9 x 10^*^ cm~^. The 
temperature of the plasma component was 0.24 (0.18- 
0.43) keV, and the abundance was 0.016 (< 0.024). Note 
that th e bulge velo city dispersion in NGC 3998 is 304 
km s~^ ((McElroy 1995), implying a virial temperature of 
~ 0.6 keV. Since NGC 3998 is an early-type galaxy and 
the inferred abundance from this fit is low, it is likely that 
this component is due to virial heating of the ISM (rather 
than starburst emission as is often observed in LLAGN). 
The power-law slope in this case was 1.84 (1.84-1.89), 
unchanged within the errors from the simple power-law 
fit. The thermal component obviously contributes to the 
spectrum most significantly at the lowest energies, with 
a 0.3-0.7 keV flux of 4 x lO'^^ erg cm'^ s-\ or ~ 13% 
of the total 0.3-0.7 keV flux. However note that in this 
band the flux estimates of the 3 EPIC detectors differ 
by ^ 8% and therefore this component should be treated 
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Table 1. Observation Parameters 



Satellite 


Bandpass (kcV) 


Date 


Exposure Time (ks) 


ASCA 


0.4-10.0 


05/10/94 


40 


BeppoSAX 


0.1-100.0 


06/29/99 


24 (LEGS), 77 (MEGS), 38 (PDS) 


XMM-Newton 


0.2-10.0 


05/09/01 


12 (MOS), 9 (FN) 




0.5 1 2 

E (keV) 



Fig. 1. — Fower-law fit to the XMM-Newton spectra. The data 
and best-fitting model are shown in the top panel and the ratio of 
data to model is shown in the bottom panel. The FN, MOSl and 
MOS2 data and residuals were marked with filled circles (black), 
filled triangles (red) and filled squares (green), respectively. 



with caution. 

In the ASCA GIS spectrum there is a small excess of 
counts near 6.4 keV which results in a "detection" of an 
Fe-K line. However the addition of a narrow Gaussian 
line (cr = 0.01 ke V) to the ASCA power -law fit only re- 
duces by 3.7 <)Terashima et al.l 12002(1 . which is only 
significant at a confidence of ~ 96% based on the f-test. 
In Table O we list the upper-limits obtained for narrow 
line emission with the energy fixed at 6.4 keV. Not sur- 
prisingly, the limit obtained when ASCA data is included 
in the fit is higher than when the XMM-Newton and Bep- 
poSAX data are fit alone. In this table we list 90% er- 
rors for one interesting parameter (A^^ = 2.7), which 
are appropriate since here there is only one parameter of 
interest (the normalization of the line), however we also 
give Ax^ = 4.6 errors for comparison with 2-parameter 
errors in the literature (or alternatively 97% 1-parameter 
errors) . We also attempted to fit a disk line model (with 
the energy fixed at 6.4 keV, inclination fixed at 0°, power- 
law emissivity slope fixed at -2 and the inner and outer 
radii fixed at 5 and 500 in Schwarzschild units), and ob- 
tained upper-limits (for Ax^ = 2.7) of 33 eV and 42 eV, 
without and with the ASCA data included, respectively. 

3.2. Optical Monitor Data 



The OM detector exposed five 1000 s images in the 
U and UVW2 filters, for a net exposure time of 5000 s 
in each band (additional exposures were taken of flank- 
ing fields which will not be discussed here). The central 
regions of the U and UVW2 images are shown in Fig- 
ure In order to properly model the surface brightness 
distributions for the two filters we also required model 
PSFs (there were no stars sufficient for this purpose in 
our observations). We initially produced model PSFs 
using SciSim. However, the current version of SciSim 
(3.0.0) does not specifically model the PSF-dependence 
of the individual filters, and also only models the PSF 
as a single Gaussian. Nevertheless, the resultant PSF 
FWHM was ~ 2.1" which is similar to that expected for 
the UVW2 filter and therefore we proceeded with this 
PSF model for that filter. In the case of the U band filter 
we instead produced a model PSF by generating double 
Gaussian images and adjusting standard deviations of 
the two Gaussians and their relative normalization until 
the encircled energy function (i.e., the integrated radial 
profile) matched the values given in the OM calibration 
file OM_PSF1DRB_0006.CCF. In the surface brightness 
fitting described below we model the nuclear point-source 
component with a Gaussian whose standard deviation 
(fj) is permitted to be a free parameter. The fitted value 
of a would then represent some combination of telescope 
jitter during the observation and also any adjustment to 
the PSF due to the inadequacies of our PSF model. 

For UVW2 data, the XMMSAS tool "omichain" pro- 
duces OM point source lists based on a 6" aperture and 
background taken from an aperture with radii 7-12" . The 
nuclear 2 120 A magnitude was 15.1 which corresponds to 
a flux density of 4.9 x 10~^^ erg cni~^ s~^A~^ (based on 
a Vega zero point and flux of 0.025 and 5.39 x 10~^ erg 
cm~^ s~^A^^^). The UVW2 images were fit with a dou- 
ble Gaussian model (with one Gaussian intended to rep- 
resent any diffuse fiux) from which it was determined that 
~ 50% of the UV flux from NGC 3998 is extended over a 
~ 6" (FWHM) scale (the individual fits typically resulted 
in Gaussian a values of ^ 1" for the nuclear component 
and ~ 3 — 4" for the diffuse component). Our spatial 
analysis implies that the flux density from within a 2" 
aperture was ~ 2.5 x 10~^^ erg cm~^ s~^A~^. Similarly, 
integrating the Gaussian flt for the nuclear component 
results in a count rate of ~ 0.5 counts s~^, which imphes 
a flux density of ~ 2.4 x 10""'^^ erg cm~^ s~^A~^ (using a 
count rate to flux conversion factor of 4.871 x 10^^^ erg 
cm~^ s~^A~^ counts""'^ s). This result is in fairly good 
agreement with the values deter mined spectrally from 
lUE data in lR.eichert et alJ l|1992D . 1.5 x 10"^^ erg cm^^ 
s~^A~^ and 1.2 x 10~^^ erg cm~^ s~^A^^ for the unre- 

^ all OM fiux conversions discussed here are based on the 
calibration given at 

[htt p: / / xmm.vilspa.esa.es/sas/documentation / watchout /uv flux.shtml 
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Table 2. Power-Law Fits to X-ray Spectra 



Detector 


Nh 


r 


-^2-10 kcV 


X^/dof 


Reference 




(1020 cm-2) 




(10~^^ erg cm~2 s^^) 






ASCA 


8.0 (6.9-9.2) 


1.89 (1.86-1.92) 


0.80 


310.9/281 


This work'^ 


BeppoSAX 


5.4 (3.5-7.8) 


1.87 (1.83-1.90) 


1.2 


222.8/227 


This work^ 


BeppoSAX, <10 keV 


5.1 (3.4-7.5) 


1.86 (1.82-1.90) 


1.1 


207.9/210 


This work 


XMM-Newton 


3.3 (3.0-3.6) 


1.88 (1.87-1.90) 


1.1 


1422.6/1338 


This work 


XMM-Newton + BeppoSAX 


3.3 (3.0-3.6) 


1.88 (1.87-1.89) 


1.1 


1652.5/1567 


This work 



Note. — Errors arc given in parenthesis for the 90% confidence interval assuming one interesting parameter (i.e., Ax^ — 2.7). 
iSee also Tcrashima et al. |2003) 
2 See also iPellegrini et al.i Ii200nlh 



Table 3. Fe-K Equivalent Width Limits 



Detector(s) EW, Ax^=2.7 EW, Ax2=4.l 

XMM-Newton < 32.7 <41.4 

BeppoSAX <38.2 <56.0 

ASCA 85.5 (11.7-159) <186 

XMM-Newton -|- BeppoSAX <25 eV < 33 eV 

XMM-Newton -\- BeppoSAX -f ASCA <33 <40 



Note. — In all fits the line energy and physical width were fixed at 6.4 kcV and 0.01 kcV (i.e., a narrow line), respectively. 




NGC 3998 XMM + BeppoSAX + ASCA 




E (keV) 



E (keV) 



Fig. 2. — Power-law fit to the XMM-Newton, BeppoSAX and ASCA spectra. The normalization of each data set was allowed to vary 
independently. The column density of a cold absorber was allowed to vary independently for each observation (but was fixed to be the same 
for each XMM-Newton detector, for example). The data are marked as follows: XMM PN, MOSl, MOS2 = filled circles, filled triangles, 
and filled squares; ASCA SIS = filled star, GIS = crosses; BeppoSAX LEGS, MEGS2-3 = circles of increasing size, PDS = unmarked (the 
only data point beyond 10 keV). The data and best-fitting model are shown in the top panel and the ratio of data to model is shown in 
the bottom panel of the left plot, and the right plot shows only the data/model ratio in the 5.0-7.5 keV region. 
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solved and extended (FWHM ~ 10") continuum compo- 
nents in the 1950-2155A region. Given the uncertainties 
(particularly in the PSF model), we note that the ac- 
tual nuclear flux density most likely lies in the range of 
2.5 — 5.0 X 10~^^ erg cm^^ s^^A~^. We therefore assume 
a flux of 3.7 X 10"^^ erg cm"^ s~^A~^ with a 50% un- 
certainty. However, this value is a factor of ^ 2 — 4 lower 
than the FOG lower limit (1.0 x IQ-^'^ erg cm^^ s-^l-i) 
determined by Fabbiano. Fassnacht. & Trinchieri (1994). 
Since the FOG UV data was overexposed in that obser- 
vation our analysis will be based on the OM flux, and we 
note that it is possible that HST observed a flare in the 
UV bandpass. 

The U band magnitude was determined from a 3" aper- 
ture (using the XMMSAS tool "omsource" which allowed 
for smaller source radii in the case of U band images). 
The U band magnitude was 14.3, which corresponds to 
a flux of 6.1 X 10~^^ erg cm~^ s~^A~^ (with the Vega 
zero-point and flux being 0.025 and 3.16 x 10"^ erg cm~^ 
s~^A~^ flux). The U band flux contains a contribution 
from the extra-nuclear emission but can of course nev- 
ertheless be considered to be an upper-limit to the nu- 
clear U-band flux of NGG 3998. We fit the central im- 
ages with a model that consisted of a Gaussian (again to 
represent a point source including jitter) and a de Vau- 
couleurs profile (oc exp — 7.669[(^)^/'* — 1.)], where Tg 
is the scale length of the profile. These fits resulted in 
re values of ~ 30 — 40" , although the precise values are 
not particularly meaningful given the uncertainty in the 
PSF and the fact that we only fit the central 30" of the 
images since we are mainly interested in the determin- 
ing the nuclear flux. The Gaussian a values were ~ 1.5", 
with the Gaussian component contributing 40-60% of the 
flux within a radius of 3" . Here the count rate integrated 
from the nuclear Gaussian component is 11 counts s~^, 
which results in a flux density of 2.1 x 10~^^ erg cm~^ 
g-i^-i_ T he U band HST FOG image was also ove rex- 
posed, and lFabbiano. Fassnacht. & Trinchieril l) 19941) re- 
port a lower-limit of 4.2 x 10~^^erg cm~^ s~^A~^ which 
is consistent with our flux. 

The current flux calibration uncertainty for OM U and 
UV magnitudes is ~ 10% (U.S. XMM-Newton Guest Ob- 
server Facility, private communication). To check the 
calibration of the U band magnitudes, we computed the 
the (total) count rate of NGG 3998 and another galaxy 
in the U band FOV, NGG 3990, and derived U band 
magnitudes of 11.92 and 13.94. These values can be 
compared to the RG3 U band magnitudes of 12.13 and 
13.81, respectively, or differences of 0.21 and 0.13. This 
implies that the flux calibration is indeed accurate to at 
least ^ 10%. As with the UV flux, the nuclear U band 
flux the uncertainty is dominated by the contribution of 
the extra- nuclear component, and therefore we assume a 
value of 3.0 x 10^^^ erg cm^^ s^^A^^ (50% of the aper- 
ture photometry value) with a ~ 30% uncertainty (i.e., 
with ItT being the offset to the model nuclear flux value 
discussed above). 

In order to supplement the OM data with additional 
optical photometric points, we derived nuclear fluxes 
from the HST WFPG2 F547M (u2uh0505t, 200 second 
exposure) and F791W (u2uh0507t, 80 second exposure) 
observations of NGG 3998. For both images the net 
count rate was computed within a 0.1" circular aper- 



ture (with the background estimated from an annulus at 
2.5-2.6"). The count rates were converted to fluxes using 
the PHOTFLAM FITS keywords (2.96 x lO-^^ erg cm'^ 
s-^l-i and 7.69 x IQ-^* erg cm'^ s'^A"^ for F791W 
and F547M, respectively). This resulted in flux densities 
of 7.4 X 10"^^ erg cm"^ s-'^A'^ and 7.9 x 10"^'^ erg cm"^ 
s~^A~^ at 7872A and 5484A, respectively. We produced 
model PSFs for these images using TinyTim and fitted 
the central 3" (i.e., a 64 x 64 pixel image) with a Gaus- 
sian plus de Vaucouleurs model. This resulted in nuclear 
flux density estimates of 6.6 x 10^^® erg cm~^ s~^A~^ 
and 7.6 x 10^^^ erg cm~^ s^^A~^, which do not differ 
by more than ~ 10% from values obtained above from 
aperture photometry, and we assume an error of 10% on 
these fluxes. 

3.3. Variability 

No statistically-significant short-term variability has 
been observed in the X-ray data, except for a marginally- 
signifi cant, ^ 12% sca.le vari ation observed in the MEGS 
data ( Pell egrini et alJl2000() . In Figure 0] we show the 
long-term 0.5-2.0 keV fight curve of NGG 3998. These 
fluxes were derived from the ASCA, XMM-Newton, and 
BeppoS AX data discussed in this paper as well as archival 
ROSAT PSPG and HRI data (in which we converted 
the count rate to the 0.5-2.0 keV flux using a power-law 
with a s lope of 1.88). The Ein stein data point from 1979 
is from iBurstein et al.l l)l997j) converted to the 0.5-2.0 
keV bandpass. We chose 0.5-2.0 keV bandpass since it 
is within t he e nergy ranges of these detectors, and, as 
shown in i)3.1l there is a negligible contribution from 
hot gas to the spectrum (as would be expected since 
NGG 3998 is not a star-forming galaxy). From this plot, 
it appears that the flux of NGG 3998 was fairly con- 
stant during the measurements but may have been a fac- 
tor of ~ 50% higher during the BeppoSAX and XMM- 
Newton observations. The Gina 2-10 keV flux for NGG 
3998 was also high, by a factor of ~ 1.9 (~ 1.3) rela- 
tive to the ASCA (BeppoSAX) fluxes, although it should 
be kept in mind that Ginga was not an imagin g tele- 
scope and had a FOV of 1° l|Awaki et al.lll99H) . Note 
that X-ray binaries are not expected to be a significant 
contributor to the X-ray luminosity of NGG 3998 since 
typical normal galaxy l uminositi es ten d to be < 10^° 
erg s~'^ l|David. Jones, fc FormanI 1199^ . in contrast to 
~ 3 X 10"'^ erg s'^for NGG 3998. More specifically, 
the B-band luminosity of NGG 3998 is ^ 1.6 x W^^eig 
s~^(based on Bt = 11.6 from RG3), which implies a 
0.5-4.5 keV l uminosity of ~ 1.6 x 10'^^ erg s~ ^ using the 
correlation in lDavid. Jones, fc Formtail 1)19921) . 

3.4. NGC 3998 X-1 

NGG 3998 X-1 lIRoberts fc Warwickll200?i | ) is an extra- 
nuclear X-ray source just outside of the RG3 (i25 ellipse of 
NGG 3998. In FigureElwe show the PN image of the cen- 
tral 4' of NGG 3998 with the position of X-1 marked. 
X-1 was detected with 441, 73, and 84 net counts by the 
PN, MOSl and M0S2 detectors, respectively, in 27.5" 
regions. Background spectra were extracted from circu- 
lar regions of the same size in a source-free region away 
from the ULX and NGG 3998. A power-law fit resulted 
in Nh = 0.5(< 5.7) x cm'^, F = 1.8(1.6 - 2.1) 
and xVdof = 25.9/32 (see Figure 0). The 0.5-2.0 keV 
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Fig. 3.— OM images of the central region of NGC 3998 in the U (left) and UVW2 (right) bands 
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Fig. 4.— Long-term 0.5-2.0 keV light curve for NGC 3998 . All 
points are plotted with errors at the level of 20% (~ 2x the varia- 
tion in flux calibration typically observed between X-ray detectors. 
See Snowden (2002). 



and 2-10 keV fluxes (averaged from the three detec- 
tors which differed from the mean by less than ~ 25%) 
were 5.1 x 10^^** erg cm^^ s^^and 8.7 x 10"-^^ erg cm^^ 
s~^. If this source is associated with NGC 3998 then 
the imphed luminosity is 2 x lO'^^erg s~^, making it an 
"ultra- luminous X-ray source" (ULX). A multi-color disk 
(MCD) black body model resulted in an unacceptable 
fit, with Nh = 0. cm-2, m„ = 0.7 keV, and xVdof 
= 53.1/32. The observed X-ray spectrum is consistent 
with ULXs, often observed with power-law spectra (e.g., 
iStrickland et al.ll200lHFoschini et alJl2002 1. However of 
course a power-law spectrum with T ~ 1.9 is also con- 
sistent with a background AGN. The XMM-New ton 2-10 
keV logN-logS shown in ' Hasinger et al.l 1)200 If) implies 
that ~ 10 background sources deg~^ would be expected 
at this flux level or higher. X-1 is 1.6' from the nucleus 
NGC 3998, and so only ~ 0.02 background sources are 
expected at this flux or higher within this area, implying 



that it is unlikely that this is a background source. Since 
we used a source-free region for the background spec- 
tra to avoid inaccurately subtracting flux due to NGC 
3998 (as would result from using a "local" background), 
however the spectra of the ULX spectra then contain 
some contamination from NGC 3998. To quantify this, 
we used our SciSim simulation of the EPIC data and 
estimate that 0.19% of the flux of NGC 3998 is contam- 
inating the X-1 source regions, or ~ 2 x 10~^^erg cm~^ 
s-i= 24% of the observed flux of X-1. 

We also show in Figure |S1 the U band image. After 
adjusting the astrometry by applying a ^ 1" shift (de- 
rived by aligning the FIRST position of the nucleus of 
NGC 3998 and two USNO-B stars), an additional USNO- 
B source, 1454-0216445, is coincident with a U band 
source and NGC 3998 X-1. The U band source is de- 
tected with a signal-to-noise ratio of 9.3 and a magnitude 
19.9 ± 0.3 (derived using the tool "omsource"; note that 
here the statistical error on the U magnitude is larger 
than the calibration uncertainties). The USNO-B 2nd 
epoch (1974.9) R and B magnitudes are 20.0 and 21.6. 
The USNO-B 1st epoch values are brighter by 0.5 and 
0.7 magnitudes, respectively. However th e uncertainty 
of US NO-B magnitudes is typically ^ 0.3 (jMonet et alJ 
12003^ . and since these magnitudes are close to the flux 
limits of the survey (V ~ 21), their uncertainties are 
probably somewhat larger, particularly in the case of the 
B band magnitude. Therefore the difference between the 
two USNO-B epochs is probably just due to statistical 
fluctuation. Similarly the unusual "concave" spectrum 
implied by the USNO-B second epoch R and B magni- 
tudes and the OM U band magnitude is largely due to 
the questionable B band magnitude. 

4. DISCUSSION 

As discussed in lPelleerini 'eral]lf2000tl . the X -ray spec- 
trum of NGC 3998 is flt weU in the 0.1-100 kcV bandpass 
with a simple power-law model that is only moderately 
absorbed {Nh ~ 3x 10^° cm^^ in contrast to the Galactic 
value of ^ 1 X 10^° cm~^). We do not find any evidence 
for large-scale (i.e., AL/L > 25%) variability on short 
(< 1 day) time scales but NGC 3998 may have varied in 
the 0.5-2.0 keV bandpass on long time scales. Our main 
goal is to assess the physics of the accretion in NGC 3998, 
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Fig. 5.— NGC 3998 PN (left) and U (right) images of NGC 3998. The PN image is marked with the RC3 ^25 eUipse for NGC 3998 
as well as the position of NGC 3998 X-1. The U band image shows the FIRST position for NGC 3998 along with USNO-B 1454-0216435 
and USNO-B 1454-0216440 after shifting by ~ 1" in order to align with the nucleus of NGC 3998 and two stars in the field. Also shown is 
the position of USNO-B 1454-0216445 which aligns with a U band source that may be associated with NGC 3998 X-1. The XMM-Newton 
X-ray position for X-1 is ll°57'"50^ 55''28™34'' J2000. The position of USNO-B 1454-0216445 is ll°57'"50.27^ -|-55''28'"34.9'' {oSset 
2.1" from the X-ray position). 



E 



NGC 3998 X- 



1 galaxies, and the 2-10 keV luminosity of NGC 3998 



1 1 —2 

erg cm 



0,5 1 2 5 

E (keV) 

Fig. 6.— Power-law fit to the PN, MOSl and MOS2 data 
(marked with filled circles, triangles, and squares, respectively) for 
NGC 3998 X-1. 



and below we compare these observed X-ray properties 
of NGC 3998 with those observed in other AGN, and 
compare the SED of NGC 3998 with Seyfert 1 and RIAF 
models. We also speculate on the constraints that our 
results can place on the geometry of the nucleus of NGC 
3998. 

4.1. Comparison with Seyferts 

The presence of broad Ha in NGC 3998 suggests 
that the nuclear region might be similar in structure 
to "typical" Seyfert 1 galaxies. The X-ray slope of 
NGC 3998 (r = 1.9) is consistent with typical Seyfert 



(2.6 X 10^1 erg s^^for i^2-io kcV = 1-1 x 10" 
s~^) is only m arginally lo w for typical Seyfert galax- 
ies (jNandra. George. Mushotzkv. Turner, fc Yaooobl 



Il997a|) . However, NGC 3998 differs in two important 
respects from typical Seyfert 1 galaxies. First, the lack 
of strong X-ray variability is in contrast with typical 
Seyferts, particularly on short time scales, since the trend 
there is for variability to increase with decreasing lumi- 
nosity i Nandra. George. Mushotzkv. Turner, fc Yaooobl 
Il997ajl . iPtak et alJ l|199aD suggested that this is tvoical 
of LINERs and LLAGN (including NGC 3998) and 
that the X-ray source regions are larger in LINERs and 
LLAGN than in typical Seyferts. One possibility is that 
the X-ray production is occurring in optically-thin flows 
(see below) rather than flares in optically-thick accretion 
disks. 

The X-ray spectrum is also remarkably featureless 
in contrast with typical Seyfert 1 galaxies, partic- 
ularly with the lack of any Fe-K emission at 6.4 
keV. The upper-limit to the equivalent width (EW) 
in the XMM-Newton and BeppoSAX joint fit is 25 
eV. Figure shows the distribution of narrow-line 
EWs r neasured in a sample of Seyfert 1 galax- 
ies l)Nandra. George. MiLshotzkv. Turner, fc Yaooobl 
ll997cD . where we also show the upper-limit for NGC 
3998. The trend in Seyfert Is is to have higher EW Fe-K 
lines with decreasing luminosity, and accordingly the lack 
of a line in NGC 3998 is conspicuous. A cautionary point 
is that Figure [7| is based on Seyfert Is with Fe-K lines 
and does not account for Seyfert Is lacking Fe-K lines. 
However, the vast majority of (X-ray bright) Seyfert 
1 galaxies appear to exhibit some level of Fc-K emis- 
sion (Nandra. George, Mushotzkv. Turner. & Yaooo^ 
Il997bl) . with very few Seyfert Is being known that lack 
a Fe-K line when the spectra are of sufBcient quality for 
Fe-K emission to be detected (T. Yaqoob, priv. comm.). 
Thus, it is likely that the accretion physics and/or 
geometry in NGC 3998 differs somewhat from typical 
Seyfert 1 galaxies. 
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Fig. 7. — Observed narrow-line Fe-K EWs as 

a function of luminosity for Scyf ert 1 galaxie s from 
iNandra. Gcorec. Mushotzkv. Turner, & Y agoobl I1997d) with 
the upper-limit for NGC 3998 also shown. 



4.2. Radiatively- Inefficient Accretion Flows 

By integrating the best-fitting power-law to the XMM- 
Newton and BeppoSAX data (T — 1.88 with normaliza- 
tion 3.8x10^"^ photons keV^^ cni^^ s^^) we calculate the 
unabsorbed 0.1-100 flux to be 6.4x 10~^^erg cm~^ s~^, or 
a luminosity of ~ 1.5 x 10^^ erg s~^. Using the bulg e ve- 



locity dispersion of 304 km s^^ from lMcElrovl l|1995|) and 
the bulge dispersion/black hole mass relation given in 
iTremainc et al. (2002) we estimate the black hole mass of 
NGC 3998 to be ~ 7x 10^ Mq. This implies an Eddington 
luminosity {Ledd = 1-25 x lO^^^) of 9 x lO^^ erg s'^, 
so L/LEdd = 2 X 10~^ (note that we get a similar result if 
we use the typical LLAGN ratio of L2_io kcV /Lboi — 0.15 
found by Ho [1999]). One possibility is that the low lumi- 
nosities in NGC 3998 are reflecting low accretion rates. 
Recent accretion modeling suggests that at low accre- 
tion rates (< a few percent of the Eddington rate; see 
Narayan, Mahadevan, & Quataert 1998) the flow will 
have low particle densities, and in turn may be optically 
thin and radiatively inefhcient. In general these models 
predict that the X-ray spectrum from an RIAF is dom- 
inated by inverse-Compton scattering (of synchrotron 
emission) and results in somewhat steep X-ray slopes (as 
opposed to the flat slopes expected from bremsstrahlung 
emission, e.g., see Mahadevan 1997). For example, in the 
case of convection-dominated accretion flows (CDAFs), 
X-ray slopes of ~ 2 are e xpected when Lx/Lsdd > 10^'' 
l|Ball. Naravan. fc Quat aert 2001). The X -ray temporal 
and spectral properties of NGC 3998 described above are 
similar to other LINERs thought to con tain RIAFs (e.g., 
iPtak et al.lll998t iQuataert et all 11999 ^ (although again 
the X-ray spectral slope is consistent with typical Seyfert 
Is also). 

4.3. Spectral Energy Distribution 

We turn now to comparison of the SED of NGC 3998 
to the predictions of various models. The data used to 
derived the SED of NGC 3998 is listed in Table H and 



is shown in Figure |H1 The X-ray data shown has been 
"unfolded" to remove the spectral response of the de- 
tectors and have been binned more so than in Figure |21 
for display purposes. The points marked with open cir- 
cles can be considered to be upper-limits since they are 
sampling radii > 10". In the near-IR the large aperture 
points are clearly dominated by galactic emission (and 
accordingly the small aperture points contain a large con- 
tribution from galactic light). However in the radio to 
mid-IR region the radiation appears to be fairly com- 
pact (i.e., unresolved to aperture sizes less than ^ 10", 
see also Knapp et al. 1996). Also variability is evident 
among several of the radio and (possibly also) mid-IR 
measurements listed in Table 0] which also implies that 
the emission is comp act. Note th at SEDs for NGC 3998 
were also presented inlFabbiano. F assnacht. & Trinchieil 
llT99l and Pellcgrii iret all l|2000j) 



however those papers 
assume the HST FOG UV flux, which is a a factor of 
^2 — 4 higher than our OM measurement. This higher 
UV flux results in a SED that resembles a typical AGN, 
while our lower UV flux results in a SED more consistent 
with the sample of low- luminosity AGN discussed in lHol 
l| 19991) . It is also evident that in the optical band, only 
the HST data with ~ 0.1" apertures are sufficient to ex- 
clude large contributions of extra- nuclear flux. Given the 
large low-mass star contribution to near-IR emission of 
early-type galaxies, this is also true for the near-IR band. 
The relative amounts of flux in the OM U band and HST 
FOG U band suggests that this contribution is at least 
90%, and likely to be much larger at longer wavelengths. 

We assessed the impact of extinction using t he ex - 
tinction equations in Cardelli. Clavt on fc MathisI l)1989h. 
iKnapD. van Driel. fc van WoerdenI 1)19851^ found that the 
upper-limit for HI absorption to the nucleus of NGC 
3998 was ~ 4 x 10^° cm~^. Including the Galactic 
contribution of ^ 1 x lO^*' cm^^ results in a conser- 
vative estimate of Nh ^ 5 x 10^° cm~^. This value 
is consistent with the X-ray observations here with the 
exception of the ASCA data, however again we note 
that the calibration of the low energy response of the 
ASCA CCD detectors is somewha t uncertain. Th is cor- 
resp onds to Ay ^ 0.22 dJe nkins fc Savagj Il974|) , and 



from ICardelli. Clavton fc Ma,this (1989) we derive that 
- 50%, 30% and 25% of the flux at UVW2, U and B, 
respectively, would be absorbed. This level of extinction 
does not impact our conclusions. 

Comparison with RIAF Models 

For accretion onto a black hole, the gravitational 
in-fall timescale can be shorter than the radiative 
timescale; the accretion energy is then stored as thermal 
energy instead of being radiated (e.g., Rees et al. 1982; 
Narayan fc Yi 1994). In the past few years, analytical 
and numerical work has shown that the physics of 
such RIAF models is far more complex than originally 
anticipated. In particular, very little mass available 
at large radii actually accretes into the black hole; 
most of it is lost to an outflow (i.e., in an adiabatic 
inflow/outflow solution, or " ADIOS") or circulates 
in c onvective motions (e.g., iBlandfor d fc BegelmanI 

| 1999t ptonc. Pri ngle. fc B egclman H 

Igumenshchev fc Abramowica 

Ouartaert fc Gruzinov^2000r iHawlev fc Balbu^ 12002 
BaUantvne'ir^!osgl2002|) . The accretion of slightly rotat 
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Table 4. Luminosity Data for NGC 3998 



V (Hz) vLi, (ergs s ^) Aperture (") 



3.2 


X 


10* 


6.5 


X 




60 


6.0 


X 


10° 


1.4 


X 


10^7 


30 


1.4 


X 


10^ 


3.3 


X 


10^7 


24 


1.4 


X 


10^ 


3.4 


X 


1037 


45 


1.4 


X 


10'' 


3.1 


X 


1037 


5 


2.7 


X 


10^ 


7.0 


X 


10^7 


6 


5.0 


X 


10^ 


9.1 


X 


10^7 


0.5 


5.0 


X 


10^ 


1.0 


X 




0.008 


5.0 


X 


10^ 


9.6 


X 


10^' 


0.007 


8.1 


X 


10^ 


1.8 


X 




6 


1.5 


X 


10^" 


2.1 


X 


10^8 


0.1 


3.0 


X 




8.6 


X 


1041 


180 


5.0 


X 


10^2 


6.6 


X 


1041 


90 


1.2 


X 


10^3 


3.7 


X 


1041 


45 


1.5 


X 


IQis 


1.0 


X 


1042 


8 


2.5 


X 


10^3 


8.4 


X 


1041 


45 


3.0 


X 


10^3 


2.1 


X 


1041 


8 


3.0 


X 


10^3 


3.8 


X 


10^1 


6 


3.0 


X 


IQis 


3.7 


X 


1041 


6 


8.6 


X 


10^3 


1.5 


X 


1042 


5 


1.4 


X 


10^'' 


5.5 


X 


1042 


3 


1.4 


X 


lO^"* 


4.0 


X 


1042 


5 


1.4 


X 


10^'' 


5.1 


X 


10^^ 


11 


1.4 


X 


10'^'' 


8.3 


X 


10^^ 


15 


1.4 


X 


10^** 


4.0 


X 


1042 


4 


1.8 


X 


10^'' 


6.8 


X 


1042 


5 


1.8 


X 


10^"' 


9.0 


X 


10^^ 


11 


1.8 


X 


10^'' 


1.4 


X 


1043 


15 


1.8 


X 


10^"' 


6.4 


X 


1042 


4 


2.5 


X 


10^'' 


9.0 


X 


1042 


3 


2.5 


X 


1014 


7.8 


X 


1042 


5 


2.5 


X 


1014 


9.9 


X 


1042 


11 


2.5 


X 


1014 


1.5 


X 


1043 


15 


2.5 


X 


1014 


7.4 


X 


1042 


4 


3.8 


X 


1014 


1.3 


X 


1041 


0.1 


5.5 


X 


1014 


1.0 


X 


1041 


0.1 


5.9 


X 


1014 


1.6 


X 


1041 


0.1 


8.7 


X 


1014 


2.5 


X 


1041 


3 


8.8 


X 


1014 


>5, 


,5 X 1040 


0.1 


1.4 


X 


1015 


1.8 


X 


1041 


2 


1.8 


X 


1015 


>4, 


,1 X 1041 


0.1 



Reference 



Rengelink et al. (1997) 
j^immelj_va^^e^^igt^^^i£ke^ ^^^4) 

Hummel (1980) 
Condon et al. (19981 
FIRST 

Condon fc Dressel (1978) 
j^imme^^a^^e^^irgt^^^i^ke^ ^^g^ 

m rnimel et a l. ( ;i982) 

Filh o. Barthel. & Ho ( 20021 

Condon fc Dressel (19781 

Hummel^ar^e^^irst^^Di^ke^ ^^8^ 

IRAS 

IRAS 

IRAS 

Willner et al. (19851 
IRAS 

Willner et al. (19851 
Willner et al. (19851 
Willner et al. (19851 
Willner et al. (19851 
Alonso-Herrero et al. (20001 
Willner et al. (19851 
^gngmo^^^^hargles ^^^2) 
Frogel et al. (19781 
2MASS Point Source Catalog 
Willner et al. (19851 
^ongmo^^^^lm^les ^^8^ 
Frogel et al. (19781 
2MASS Point Source Catalog 
Alonso-Herr ero et al. (20001 
Jo iner et aL (19851 
^ongmo^^^^ha^les ^^8^ 
Frogel et al. (19781 
2MASS Point Source Catalog 
This work 
This work 

^^bbiano^^|Ssn^|Ch^^^^inchifin ^^9^ 
This work 

^^bbiMio^^|Ssnacht^^^Wnchign ^^94) 
This work 

^^bbimio^^ssnach^^^^inchign ^^94) 



ing gas also may result in a substantially reduced accre- 
tion rate relative to the Bondi rate LProea fc BegelmanI 
l200.1ft . Following a proposal due to Blandford & Begel- 
man (1999), we can parameterize the radial variation of 
the gas density in the flow with a parameter p, where 
p OC 7--3/2-I-P (equivalently, the accoretion rate varies 
with radius as m oc r^); p « 1/2 — 1, rather than 
p — Q (as in the original analytical advection-dominated 
accretion flow [ADAF] models of Narayan & Yi 1994), 
is favored by the simulations, implying much lower gas 
densities close to the black hole. 

Figure |Ht s hows several RIAF models for the SED 
of NGC 3998. iQuataert fc NaravanI l|1999n showed that 
there are significant degeneracies in comparing theoreti- 
cal RIAF spectra to observed data so these results should 
be interpreted as representative models, not unique 
"fits." The solid line in|SjL is an ADAF with an accretion 
rate of 3 x 10~^ (hereafter accretion rates are given in Ed- 
dington units) and 5 = 0.01 (= fraction of the accretion 
energy transfered to electrons). The dotted line is the 
ADAF model with an outer thin disk starting at a transi- 



tion radius rt = 300i?s {Rs = the Schwarzschild radius) . 
The dot-dashed line is an ADIOS model with 6 — 0.3 
and p — 0.4 and a similar accretion rate at the outer 
edge of the flow, also with an outer thin disk (the so- 
called multi-color disk (MCD) model) with rt — SOOi?^. 
Note that ADIOS models result in similar spectral mod- 
els to CDAF models. The thin disk (for both models) has 
an outer radius of lO^i?^, an inclination of 60 degrees, 
and TO w 3 X 10"'^ (hereafter all to values are quoted in 
Eddington units). Given theoretical uncertainties about 
electron heating, values of(5<<ltoi5~l are plausible 
(e.g., Quataert fc Gruzinov 1999), and the two models 
chosen in fact span this range. These models flt the IR- 
X-ray emission reasonably well, except of course for the 
large contribution of stellar flux at near-IR wavelengths. 
The (nuclear) IR emission is dominated by the (putative) 
thin disk component while the UV- X-rays are dominated 
by the RIAF, and in fact rt and the accretion rate at rt 
were chosen to approximately match the IR flux. The 
smaller rt value could have been chosen which would have 
resulted in a larger thin-disk flux, and correspondingly a 
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smaller m would have been chosen. Therefore the values 
of m and rt are only constrained to within ~ an order of 
magnitude. The X-ray emission in these RIAF models 
is dominated by inverse-Compton scattering which re- 
sults in approximately in a power-law continuum. The 
ADIOS model has a lower m and density near the black 
hole, but a larger temperature. In that situation, the 
Compton peaks are more 'separated' (since the energy 
gain per scattering is proportional to temperature) and 
the peaks are more distinct. This is a generic difference 
between high m ADAF type models and lower m models. 
Since the X-ray spectrum of NGC 3998 is fairly feature- 
less lower temperature models are more consistent with 
the X-ray data. Note that in general different thermal 
line emission characteristics are predicted for these dif- 
ferent type s of models, but these lines are not detectable 
at this TO ijPerna. Ravmond. fc Naravanl [20G0). Neither 
model accounts for the radio emission, in spite of the 
fairly large outer disk radius of lO^i?^. The most restric- 
tive constraint in this regard is the VLBI 5 Ghz point 
which has a beam size of only 7 mas, or ~ 0.7 pc, and is 
clearly nuclear in origin. 

AGN Spectral Models 

In Figure IHt) we also plot the median SEDs for 
radio- load and radio-quiet AGN from Elvis et al. ( 1994), 
rescaled to the flux of NGC 3998 at = lO^* Hz. The 
QSO SEDs grossly overestimate the blue-UV emission, 
i.e., NGC 3998 is lacking the "big blue bump" that 
is considered to be due to an optically-thick accretion 
disk. However, this is not unexpected since the accre- 
tion rate is much lower than that of a typical QSO, and 
MCD temperatures scale as rh^^^ l|Frank. King. Sz Raind 
I2002tl . We also show an approximate Comptonized thin- 
disk model that fits the X-ray data from NGC 3998 
without over-predicting optical or IR flux. This ther- 
mal Comptonization model (implemented in XSPEC as 
"compTT" ; see Titarchuk 1994, Hua & Titarchuk 1995) 
uses the Wien approximation to black body flux for com- 
putational efficiency, however this should be a reasonable 
approximation. In this model there are 5 parameters: 
the overall normalization, the temperature of the Wien 
component, the temperature and optical depth (r) of the 
hot plasma, and a geometry switch to select a spherical 
or disk geometry. In general inverse-Compton spectra 
tend to be power-laws with a slope of — where A is 
the mean amplification of the p hoton energy per scatter- 
ing ( Rvbicki fc Lightmanlll979l) . Therefore it is not sur- 
prising that for a given Wien and plasma temperature 
combination, r could be adjusted to result in our ob- 
served power-law slope and hence give a good fit. The fits 
generally resulted in a power-law extrapolated from the 
peak of the Wien component, and we therefore initially 
fixed the temperature of the Wien spectrum at 10^"^ keV 
(~ 10'^ K) to see if the power-law could be extended to the 
optical bandpass. W e tried hot plasma temp eratures of 
100 and 200 keV Csee lPietrini fc Krolikl[l995() . and fitted 
for the model normalization and r as free parameters. 
The best-fit optical depths were 0.3 and 0.1 for a disk 
geometry and plasma temperatures of 100 and 200 keV, 
respectively. With a spherical geometry the best-fit op- 
tical depths were 0.5 and 1.0 for plasma temperatures of 
100 and 200 keV. Clearly the value of r is not meaningful 
since our data do not constrain the geometry or tempera- 



ture of the hot plasma. These fits (simultaneously to the 
XMM-Newton and BeppoSAX data) resulted in x^/dof 
= 1653-1654/1567, or basically equivalent to the pure 
power-law fit. The fit with a plasma temperature of 200 
keV and a disk geometry is shown in Figure |H1 We com- 
puted a MCD spectrum (with inner radius = 3i?s and 
outer radius = 10^i?s) consistent with the Wien compo- 
nent by assuming an inclination angle of 60 degrees and 
adjusting to, which resulted in m = 4.1 x 10~^. Such a 
model accounts for the spectrum of NGC 3998 from the 
optical to ~ 100 keV, i.e., the optical flux is consistent 
with the extrapolation of the power-law observed in the 
X-ray band, particularly if we allow for variability in the 
FOG F480LP flux (the offset between the PDS point in 
Figure |S1 and the Comptonized spectrum model is due 
in part to the fact that the X-ray data were unfolded 
assuming the simple power-law spectrum extending to 
100 keV). However this model requires that a negligible 
fraction of the emission longward of the optical is due to 
the AGN in NGC 3998. This may be unhkely given the 
compactness of the radio and mid-IR emission. 

Of course there is some degeneracy between the incli- 
nation angle and accretion rate of the disk in this sce- 
nario, with an angle of 0° and m ~ 2.0 x 10^^ also 
being consistent with the Wien component. However 
the luminosity of the MCD disk component is not ar- 
bitrary. The MCD spectrum is given by the sum of 
blackbody spectra with temperatures varying as r~^^^ 
where r is radius of each poi nt on the surface of the disk 
i Frank. King, fc Ra inc 2002'). This results in a spectrum 
similar to a single blackbody emitting at the peak tem- 
perature of the disk, and for given black hole mass the 
peak temperature scales as to^^"*. Of course the lumi- 
nosity of a blackbody emitting over a given surface area 
is proportional to i/'^, where hiy ^ 3kT. Therefore we 
could not match the luminosity of the MCD component 
for other disk temperatures by simply adjusting the ac- 
cretion rate since the expected luminosity of the Wien 
component would move in the log Ljj — v plane along a 
slope of 4 rather than 0.9. This makes the Comptonized 
MCD model inconsistent with the observed X-ray power- 
law if the peak disk temperature is below ~ 10~^ keV 
(i.e., the Wien component normalization would be larger 
than the MCD expectation), and also at higher temper- 
atures this model is only consistent with the data if the 
MCD component is attenuated (e.g., by inclination ef- 
fects). Also note that the optical (UV) flux would not 
be due to the accretion disk for disk temperatures in 
excess of ^ 10~^(10~'^) keV. Conversely, assuming that 
the optical and UV data are due to a thin disk requires 
that TO < 10"^. Finally, note that a standard a disk is 
both stable and optically thick when m is in the range 
discussed here (10~^ < m < 10~^). 

Jet Models 

In addition to a very low accretion rate, another pos- 
sibility is that both the radio and X-ray flux from NGC 
3998 are due to a jet, as has been sugges ted recently for 
LLAGN such as NGC 4258 l|Yuan et al.l [20021. I C 1459 
llFabbiano et"ani200l and NGC 4594 llPellegrini'ei~ai] 
2003). Along these lines we show a simple broken power- 
law model in Figure |H1 where we found that the ra- 
dio data were consistent with a slope of ~ (i.e., with 
vL^ oc z^^) and break frequency of ^ 10^^ Hz. This is 
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1q8 ^qIO iq12 ^qU ^qIS ^ q20 ^ q8 ^ qIO ^q12 , q14 ^q16 ^qIS ^ q20 

V (Hz) 1/ (Hz) 

Fig. 8. — The SED of NGC 3998 plotted with (left panel) two RIAF models and (right panel) the radio-load (red solid line) and 
radio-quiet (green dotted line) mean QSO SEDs from lElvis etaT] tWM> . The U and UV points (at u ~ 10^5) are due to the OM, and 
the X-ray points are due to the spectral fits discussed in the text (note that the data points have been "unfolded" to account for the 
spectral responses of the detectors and hence are somewhat model dependent). The turn-over in the X-ray data at low energies is due to 
absorption (with Nh ~ 3 x 10^" cm~^). The X-ray data points with lower luminosities are due to ASCA (see text). The RIAF models are 
computed with an inner ADAF (green dotted line; m = 3 X 10~^ in Eddington units, <5 = 0.01) and an inner ADIOS (blue dot-dashed line; 
m = 3 X 10~^,<5 = 0.3, p = 0.4), both with an outer thin disk with a transition radius of 300 Rs and outer radius of lO^Rs- We also plot 
the ADAF model with no outer thin disk for comparison (red solid line). The right panel also shows a Comptonized disk spectrum (black 
bold dotted line). The Comptonization model used is based on a Wien spectrum approximation, and we show the equivalent multi-color 
disk model for comparison (blue dot-dashed line; the parameters were m = 4.1 X 10~®, and inclination angle = 60°). A simple broken 
power-law model with a radio slope of 0.0 and an IR-to-Xray slope of 0.9 is also shown in the right panel as a representative synchrotron 
model. 



obviously an overly simplistic model but here the radio 
slope is due to optically-thick synchrotron emission and 
the X-ray emission is due to optically-thin synchrotron 
emission, modeled by IMarkofF. Falcke. &: Fended lj200lD 
as being due to the outer, "post-shock" portion of the 
jet. The Markoff-Falcke model also predicts a mid-IR 
"bump" due to synchrotron emission from the inner, 
"pre-shock" jet which may be consistent with the mid- 
IR fluxes observed in NGC 3998. Alternatively, the syn- 
chrotron emission may break to a flatter slope and the 
X-ray is again also due to inverse-Compton emission, 
and the mid-IR and/or optical-UV flux is due to a stan- 
dard a disk IpVIarkoff ct al. 2003). Clearly these models 
would also be consistent with the SED of NGC 3998 and 
may imply accretion rates possibly orders-of-magnitude 
higher than the RIAF and Comptonized a-disk models 
since a large fraction of the accretion energy would be 
released kinetically in the jet. However, as was the case 
in the SED modeling discussed above, any a disk in this 
scenario would likely still need to be truncated to avoid 
over-predicting the optical and UV fluxes. Also, the lack 
of any short-term, large-amplitude X-ray variability may 
argue against the X-ray emission being dominated by a 
(relativistic) jet. 

4.4. Constraints on the Nuclear Geometry of NGC 3998 

The Fe-K lines in Seyfert 1 galaxies are assumed to 
originate from reflection from a thin disk and/or mat- 
ter further out such as the putat ive "torus" in the stan- 
dard model ijYaaoob fc Padmaiiabhan.2002 ). If the Fe-K 
emission is dominated by a thin disk, then the maxi- 



mal solid angle as seen from the X-ray source would be 
^ 27r, and the solid angle of a toroidal geometry of matter 
most likely would not exceed this. If we assume that the 
EW ~ 100 eV Fe-K line typically observed in Seyfert 1 
galaxies is produced by a reflection from optically-thick 
material with Q, 2tt, then our upper-limit of 25 
eV implies tha t »/27t ~ 0.2 5 in NGC 3998. Interest- 
ingly, Pellegrin i et al.l (|200nr i report a similar result of 
n/2Tr < 0.21 for either a neutral or ionized Compton re- 
flection hump based on reflection model fits to the Bep- 
poSAX spectrum of NGC 3998. We also computed an 
upper-limit for a disk line model at 6.4 keV and obtained 
- 30 eV. A L2-10 keV = 10''^ erg s'^ Seyfert 1 would 
typically exhibit a EW = 300 eV disk line, which im- 
plies Jl/27r < 0.1. This suggests that the solid angle of 
any optically-thick material in NGC 3998 is less than at 
most 25% of that typically observed in Seyfert 1 galaxies. 

We can roughly constrain the size of an outer thin disk 
using the line information. Assuming a disk geometry 
with the X-ray source at a height h above the center 
of the disk with inner and outer radii and rout, the 
solid angle subtended by the disk is given by = 27r[(l-|- 
(I^)2^)-l/2_(^^(r^)2)-l/2]^ The X-rays are produced 
in the inner part of the RIAF, implying h ~ 10 — 30i?s. It 
is therefore likely that rt >> h. This implies that ~ 
h{-^ ^) and therefore rt > ^, or > 100 - SOORs 

m out 

for 1^ < 0.1. This is consistent with Figure |H1 where 
we assumed ri„ = rt 300 Rs- A possible caveat to this 
analysis is that NGC 3998 may contain an optically-thick 
disk with ri„ — 3Rs but whose surface layers are fully 
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ionized. This could result in little or no Fe-K em ission 
^Done Navakshinll2f)f)TllR^,]la,ritvrie R,nsai2nnl. 

4.5. The Nature of the ULX NGC 3998 X-1 

If this is source is associated with NGC 3998, then the 
optical luminosity of the source (M^j ~ -11) implies sev- 
eral to tens of supergiant stars or a cluster of hundreds 
to thousands of low-mass stars, if the optical flux is stel- 
lar. If the optical variability observed between the two 
USNO-B epochs (separated by of order 25 years) were 
real it would argue against this, although again the vari- 
ability is probably not statistically significant. In addi- 
tion, the high U-band absolute magnitude [Mu = —10.9) 
and very blue color imply that any such cluster would be 
composed of early-type stars. This rules out a globular 
cluster with late-type stars which would be expected in a 
early type galaxy. Using logfy/F„„, = loeffn .=^-■^ w.v) 
+ R/2.5 + 5.7 ('Honschemei er et al.ll2nn3l: iSzokolv et all 
|2004; Norman et al. 2004), we derive \og Fx I Fopt = 0.4 
fo r the ULX which is typical of AGN found in the CDF- 
S l|Szokolv et alJl2(fflllNOTman et alJl200l and is much 
larger than values typical of stars (logFx/Fopt « —1) 
and somewhat smaller than values typical of low-mass 
X-ray binaries {log Fx /Fopt > 1; Bradt & McClintock 
1983; although realistically we would be observing an X- 
ray binary in a cluster). The most likely identification 
for this source is therefore a background AGN, although 
as stated in i l3.4l the detection of a background source 
at this flux was unlikely.. 

5. SUMMARY 

We have determined the tightest constraint to date on 
Fe-K emission in NGC 3998, specifically EW < 25 eV for 
a narrow line at 6.4 keV. Our analysis of XMM-Newton 
data combined with archival ASCA and BeppoSAX data 
shows that the spectrum is fit well with a simple power- 
law with r = 1.9 and a small amount of intrinsic absorp- 
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tion. The XMM-Newton OM UV flux measurement is 
consistent with an extrapolation of the X-ray power-law, 
suggesting that the UV to X-ray spectrum of NGC 3998 
is consistent with a power-law over 5 decades in energy. 
The SED of NGC 3998 is consistent with both RIAF 
models and a Comptonized thin-disk model (i.e., as ex- 
pected to be found in Seyfert galaxies), both of which 
can account for the UV to X-ray spectrum of NGC 3998. 
None of the accretion models account for the compact ra- 
dio emission in NGC 3998 which is probably dominated 
by a jet. The "pure" thin-disk model also does not ac- 
count the mid-IR emission, although a truncated thin- 
disk surrounding an inner RIAF fits the mid-IR to X-ray 
spectrum well (with the near-IR to optical region dom- 
inated by extra- nuclear flux). If the accretion is domi- 
nated by an RIAF then the transition radius is > 300i?s 
and the m is ~ 3 x 10-3 (at ~ 300i?s in the models with 
p > 0). In either case the inner radius of the thin disk 
is probably > 100i?s in order to account for the lack of 
X-ray reflection features. If the accretion is dominated 
by a Comptonized thin disk that extends to a radius of 
3i?s, then the accretion rate in the disk is in the range 
10"^ — IQ-^ in Eddington units and the lack of Fe-K 
emission may be due to an ionized accretion disk. Fi- 
nally, a jet model may be consistent with the radio and 
X-ray flux, and jet emission may also be responsible for 
the mid-IR and optical/UV flux. 
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